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List  of  Symbols 

angular  momentum 
total  pressure  coefficient 
static  pressure 
radius  of  curvature 
radial  position 
axial  velocity 
tangential  velocity 
meridional  velocity 
free  stream  velocity 
outlet  flow  angle 
section  turning  angle 
flow  coefficient 
propeller  angular  velocity 


Subscripts  and  Superscripts 

Ap  area  of  propeller 

E  denote  propeller  exit  station 

I  denote  propeller  inlet  station 

r^  radius  of  inner  streamline 

r  radius  of  bounding  streamline 

o 

r 

P 


1 

i 


radius  of  propeller 


the  propeller  designer  establishes  certain  performance  objectives  such 
as  the  required  thrust,  torque,  vehicle  speed,  or  other  so-called  "design" 
conditions.  With  this  information,  the  flow  field  is  calculated  and  the  proper 
propeller  geometry  is  determined.  On  the  other  hand,  the  direct  method 
deals  with  the  analysis  of  a  propeller  with  a  known  geometry.  For  this 
case,  the  flow  field  and  propeller  performance  data  are  calculated  for 
design  and  off-design  operating  points. 

The  open  propeller  is  unique  in  comparison  to  other  types  of  turbo¬ 
machines.  As  the  name  indicates,  this  is  an  "open”  machine,  i.e.,  the 
flow  field  lacks  a  well-defined  boundary,  such  as  the  walls  of  a  turbine 
or  compressor.  However,  it  is  possible  to  consider  the  flow  field  for 
the  open  propeller  as  being  the  region  constrained  within  a  potential 
flow  streamtube.  The  bounding  streamsurface  of  this  streamtube  is 
actually  located  infinitely  far  from  the  propeller;  however,  it  is 
assumed  that  there  is  a  finite  distance  at  which  the  flow  field  is 
relatively  unaffected  by  the  presence  of  the  propeller.  This  distance 
is  assumed  to  be  at  least  eight  times  the  propeller  radius.  The  bounding 
streamline  at  this  distance  serves  as  a  flow  boundary  in  this  analysis. 

In  this  manner,  the  streamtube  becomes  analogous  to  the  compressor  or 
turbine  duct. 

A  method  which  is  commonly  applied  to  calculate  the  flow  field  in 

an  axisymmetric  turbomachine  is  the  streamline  curvature  method  (SLC) . 

This  method  has  been  used  for  some  time  for  the  design  and  analysis  of 

1  2 

various  types  of  turbomachines.  As  described  by  Novak  ,  Frost  ,  and 
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Davis  ,  this  numerical  technique  has  been  successfully  used  for  the 

A 

indirect  analyses  of  compressors  and  turbines.  Also,  McBride  has 
applied  this  streamline  curvature  method  to  the  Indirect  problem  of  an 
open  hydrodynamic  propeller.  However,  very  little  work  has  been  done 
using  the  streamline  curvature  method  for  the  direct  open  propeller 
problem. 

It  is  the  objective  of  this  paper  to  discuss  a  method  which  can  be 
used  to  solve  the  direct  open  propeller  problem.  The  streamline  curvature 
technique  used  to  establish  the  streamlines  in  the  propeller  flow  field 
is  discussed  in  detail. 


II.  Streamline  Curvature  Method 

Basic  Equations 

The  major  equations  used  in  the  streamline  curvature  method  of 
analysis  are  derived  from  the  principles  of  conservation  of  mass,  momentum, 
and  total  energy.  The  equations  which  are  used  are  generalized  for  an 
axisymmetric  flow  field  and  allow  for  curvature  of  the  streamlines.  The 
working  fluid  in  this  analysis  is  assumed  to  be  incompressible,  inviscid 
and  steady.  Thus,  the  equations  are  derived  with  these  assumptions  in 
mind. 

One  of  the  fundamental  equations  used  in  the  analysis  is  obtained 
from  the  principle  of  conservation  of  mass  within  a  streamtube.  The 
continuity  equation  states  that  for  a  steady  flow, 

pV  dA  =  0  .  (1) 

c.s. 


This  integral  is  evaluated  over  a  control  surface  where  V  is  the  velocity 


-5- 


18  August  1977 
HJP : j  ep 


component  normal  to  the  element  of  surface  area  (dA) .  By  using  a  stream- 
tube  as  the  control  volume,  there  Is  no  transport  of  mass  normal  to  the 
streamsurface.  Thus,  the  continuity  equation  for  an  incompressible  flow 
within  a  streamtube  can  be  simplified  to 

•out 

2tt  V  r  dr  =  0  ,  (2) 

in  A 

where  the  integral  is  evaluated  from  some  inlet  station  to  an  exit  station 
along  the  streamtube.  Because  the  areas  represented  by  this  form  of  the 
continuity  equation  are  located  in  radial  planes,  the  component  of 
velocity  is  necessarily  in  the  axial  direction  (V  ) . 

For  this  axisymmetric  analysis,  it  is  necessary  to  use  only  the 
components  of  total  velocity  in  two  directions.  The  component  of  velocity 
which  is  tangent  to  the  streamline  and  is  projected  onto  the  meridional 
plane  is  called  the  meridional  velocity  (Vffi) .  This  component  of  total 
velocity  is  related  to  the  axial  component  by  the  cosine  of  the  streamline 
angle  (a) .  The  second  component  of  total  velocity  is  known  as  the 
tangential  velocity  (Vg),  and  it  is  located  in  the  circumferential 
direction,  which  is  perpendicular  to  the  meridional  plane.  These  velocity 
components  as  well  as  the  computational  coordinate  system  are  illustrated 
in  Figures  la  and  lb. 

Meridional  curvature  of  the  streamlines  in  the  flow  tends  to  disrupt 
the  equilibrium  along  a  streamline.  To  satisfy  the  principle  of 
conservation  of  momentum,  a  radial  static  pressure  gradient  develops. 

There  are  really  three  contributing  factors  in  the  development  of  this 
pressure  gradient.  The  first  term  is  indeed  due  to  the  moving  particle 
being  subjected  to  the  streamline  curvature  in  the  meridional  plane. 
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The  magnitude  of  this  term  depends  greatly  on  the  radius  of  curvature 

of  the  streamline  at  a  specific  location  (R  ).  The  second  term  is 

c 

directly  related  to  the  centrifugal  force  which  the  rotating  fluid 
experiences  upon  passing  through  the  propeller  blades.  The  final  term 
is  proportional  to  the  convective  acceleration  of  the  fluid  particle 
as  the  flow  area  either  converges  or  diverges.  The  combination  of 
these  three  terms  is  known  as  the  radial  equilibrium  equation  and  is 
expressed  as 


dP 

P 


V2  V2 

m  ,  ,  e  . 

—  dn  +  —  dr 
R  r 


DV 

Vm  DU-  ds 


(3) 


It  can  be  seen  that  for  an  axial  flow  turbomachine  with  no  streamline 
curvature  and  no  variation  in  flow  area.  Equation  (3)  reduces  to  the 
simplified  version  of  the  radial  equilibrium  equation. 

Because  total  energy  must  be  conserved  along  a  streamline,  the  total 
energy  between  two  points  on  a  streamline,  upstream  or  downstream  of  the 
propeller  but  not  across  the  propeller,  can  be  written  in  terms  of  the 
total  pressures 


?1  +  l/2p  V2  =  Pn  +  l/2p  V2  ,  (4) 

where  the  subscript  (1)  indicates  a  station  of  known  total  pressure  and 
the  subscript  (n)  indicates  an  arbitrary  station.  This  equation  is  seen 
to  be  Bernoulli's  equation  along  a  streamline  for  an  inviscid  flow.  A 
simplification  to  this  equation  is  possible  by  expressing  the  total 
velocity  (V^)  in  its  components  and  by  realizing  that 
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where  Pp^  is  the  absolute  static  pressure  on  the  innermost  streamline 
for  a  given  station  in  the  flow  field.  The  integral  of  the  radial 
equilibrium  equation  yields  the  static  pressure  difference  between  the 
inner  streamline  (r)^)  and  an  arbitrary  outer  streamline  (r)Q) .  Using 
Equation  (5)  and  rearranging  Equation  (4),  an  expression  for  the  meridional 
velocity  is  obtained: 

rn 

l/2p  V *  =  P  +  l/2p  V*  -  Pn.  -  dP  -  l/2p  V*  .  (6) 

N 

Normally,  this  form  of  the  energy  equation  cannot  be  solved  since  there  are 
two  unknowns,  namely,  Pp^  and  Vm-  However,  by  using  continuity  and 
eliminating  the  meridional  velocity  term,  the  static  pressure  on  the 
inner  streamline  (Pp^  can  be  computed.  Once  this  value  is  known,  a 
new  meridional  velocity  profile  which  satisfies  the  energy  equation  can 
be  calculated  for  each  station  in  the  flow  field  using  Equation  (6). 

The  three  equations,  continuity,  radial  equilibrium,  and  energy,  are 
the  basis  for  the  flow  field  solution  obtained  by  the  streamline  curvature 
method . 

Boundary  Conditions 

The  streamline  curvature  method  uses  an  iterative  procedure  to  obtain 
a  flow  field  solution  which  simultaneously  satisfies  the  three  basic 
equations.  The  iteration  procedure  can  be  long  and  laborious.  For  this 
reason,  an  indirect  streamline  curvature  computer  program  developed  by 
McBride  ^  has  been  modified  to  solve  the  direct  problem  of  the  open 
propeller.  A  flow  diagram  of  the  major  computational  steps  is  shown 


in  Figure  2. 
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As  the  first  two  blocks  in  the  flow  diagram  indicate,  there  are 
two  initial  boundary  conditions  which  must  be  known  before  the 
computational  process  can  proceed.  The  first  boundary  condition  is 
that  the  flow  energy  must  be  known  at  some  reference  station  far 
upstream  of  the  propeller.  This  flow  energy  is  usually  specified  in 
terms  of  the  components  of  total  velocity.  Knowing  the  meridional 
and  tangential  velocity  distributions  along  some  reference  station  line, 
the  total  energy  along  any  streamline  can  be  calculated  and  must  remain 
constant  until  it  is  changed  by  the  propeller.  This  change  in  energy 
through  the  propeller  is  the  other  boundary  condition.  In  the  direct 
analysis,  the  energy  along  a  streamline  coming  out  of  the  propeller 
is  related  to  the  angle  at  which  the  flow  exits  from  the  propeller  blade 
at  each  streamline  position.  These  exit  flow  angles  must  be  known 
before  any  estimation  of  the  energy  change  through  the  propeller  is 
made. 

The  determination  of  these  exit  flow  angles  is  the  heart  of  this 
computational  method.  The  ideal  outlet  flow  angle  is  simply  the  fixed 
angle  between  the  tangent  to  the  camberline  at  the  trailing  edge  of 
a  blade  section  and  the  exit  axial  velocity  vector.  This  angle  is 
shown  in  Figure  3.  This  ideal  outlet  flow  angle  however,  is  not  the 
angle  at  which  a  real  fluid  exits  from  the  blade  row. 

The  calculation  of  the  deviation  term  which  allows  for  real  flow 
effects  is  a  necessary  boundary  condition  for  the  direct  problem.  A 
method  for  estimating  a  deviation  angle  will  be  discussed  in  a  later 
section.  For  the  present,  it  is  assumed  that  the  measured  outlet  flow 
angles  and  the  correct  deviation  terms  have  been  computed  and  are  used 
as  the  second  boundary  condition  for  the  analysis. 


-9- 


18  August  1977 
MJP :  j  ep 


Computational  Procedure 

As  the  flow  diagram  In  Figure  2  indicates,  the  initial  velocity 
distributions  at  the  reference  station  must  be  transferred  to  the 
downstream  flow  stations  so  that  an  initial  flew  field  can  be  established. 
This  transfer  is  carried  out  in  three  steps.  First,  the  meridional 
velocity  profile  is  transferred  throughout  the  flow  field  by  means  of 
the  continuity  equation.  Second,  the  reference  tangential  velocity 
distribution  is  transferred  through  the  flow  field  to  the  propeller 
inlet  station  by  means  of  conservation  of  angular  momentum.  Since  the 
propeller  changes  the  angular  momentum  along  a  streamline,  the  tangential 
velocity  profile  at  the  propeller  exit  station  must  be  obtained  to  allow 
for  further  calculations.  This  is  where  the  second  boundary  condition 
enters  the  analysis.  From  Figure  3,  it  can  be  seen  that  the  outlet  flow 
angle  (62)  is  related  to  the  tangential  velocity  by  the  equation 

Vq  =  wr  -  VAz  tan  $2  ,  (7) 

where  wr  is  the  angular  velocity  of  the  propeller  and  is  the  axial 
velocity  at  the  propeller  exit  station.  Knowing  the  outlet  flow 
angles  enables  the  values  of  tangential  velocity  to  be  computed  at  the 
propeller  exit  station.  Once  the  angular  momentum  at  this  station 
is  calculated,  the  principle  of  constant  angular  momentum  allows  for  the 
third  step  to  be  carried  out;  the  transfer  of  tangential  velocity 
downstream  of  the  propeller. 

After  the  initial  meridional  profile  has  been  transferred  downstream, 
the  initial  streamline  locations  can  be  determined.  Locating  the 
streamlines  is  accomplished  by  using  a  streamline  spacing  function. 
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This  function  assigns  a  specific  percentage  of  the  total  mass  flow  to  each 
streamline,  thereby  establishing  a  unique  radial  position  for  each  stream¬ 
line  along  every  station  in  the  flow  field.  These  points  are  then  fit 
with  a  line  which  represents  the  individual  streamline. 

At  this  point,  the  flow  field  has  been  established  based  on  the  two 
boundary  conditions  and  continuity.  The  solution,  however,  must  also 
satisfy  radial  equilibrium  and  total  energy  requirements.  The  pressure 
gradient  given  by  the  radial  equilibrium  equation  (Equation  (3))  is 
computed  at  every  station  from  the  inner  streamline  to  the  bounding 
streamsurface.  To  obtain  the  static  pressure  on  the  inner  streamline, 
it  is  necessary  to  combine  the  energy  equation  (Equation  (6))  and 
continuity  equation  (Equation  (2)).  The  stations  upstream  of  the 
propeller  are  calculated  based  on  the  known  total  pressures  at  the 
reference  station.  The  static  pressure  at  stations  downstream  of 
the  propeller  require  the  calculation  of  total  pressures  (static 
plus  dynamic)  along  a  streamline  at  the  propeller  exit  station.  This 
calculation  is  obtained  from  an  angular  momentum  change.  Angular 
momentum  is  defined  along  a  streamline  as  the  tangential  component 
of  velocity  multiplied  by  its  radial  distance.  Knowing  the  angular 
momentum  before  and  after  the  propeller,  the  total  pressure  in  coefficient 
form  at  the  exit  station  is  given  by. 


Rtt  ^ 

V  ■  cpt  +  rrr  (ame  -  "V  (8 

E  I  r  00  p 

where  Cpg  is  the  total  pressure  coefficient  at  the  propeller  exit  station. 
This  equation  yields  the  increase  in  total  pressure  across  the  propeller 


along  a  streamline. 
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Having  computed  the  total  pressure  at  the  exit  station,  it  is  now 
possible  to  calculate  a  new  meridional  velocity  profile  for  every  down¬ 
stream  station  which  satisfies  total  energy  conservation  (Equation  (6)). 
This  new  velocity  profile  is  then  integrated  for  a  new  mass  flow 
distribution.  Using  the  streamline  spacing  function  and  this  new 
distribution  of  mass  flow,  new  streamline  locations  are  calculated  and 
the  entire  set  of  streamlines  is  repositioned. 

The  convergence  criterion  for  this  iteration  procedure  is  related 
to  the  amount  of  shifting  that  any  particular  streamline  experiences.  If 
the  shift  of  every  streamline  is  within  prescribed  limits,  the  solution 
is  said  to  be  converged  and  the  iteration  cycle  is  complete.  If  the 
shift  of  any  streamline  is  not  within  the  specified  tolerance,  the 
procedure  continues  its  iteration  cycle  as  shown  by  the  flow  diagram  in 
Figure  2.  Upon  convergence,  the  flow  field  simultaneously  satisfies 
continuity,  conservation  of  momentum,  and  conservation  of  total  energy. 

Once  the  flow  field  has  been  calculated,  the  performance  parameters, 
such  as  thrust  and  torque,  can  be  computed  based  on  energy  and  momentum 
considerations.  The  torque  coefficient  is  defined  as 


fr0  2 7i  AMg  vmE  r  dn  frQ  ztt  aMj.  vm]; 

=  j  1/2  Vi  A  r  1/2  A 

'  r.  00  p  p  'r.  00 


r  2tt  AMt  VmT  r  dn 


This  torque  coefficient  gives  an  indication  as  to  the  size  requirements 
for  the  powerplant  to  be  used  in  conjunction  with  a  specific  propeller. 

A  measure  of  the  forward  force  produced  by  the  propeller  is  in  the 
form  of  a  thrust  coefficient.  The  thrust  coefficient  can  be  obtained 
by  means  of  actuator  disc  theory.  In  this  analysis,  the  rotating  blade 
row  is  modeled  by  a  solid  actuator  disc  and  the  momentum  equation  yields 
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the  resultant  force  in  coefficient  form  as 


r  dr 

l/2p  V‘  A 

00  p 


(10) 


III.  Outlet  Flow  Angles 
Calculation  of  Deviation  Angles 

The  propeller  outlet  flow  angle  is  defined  as  the  angle  between  the 
tangent  to  the  blade  section  camberline  at  the  trailing  edge  and  the  exit 
axial  velocity  vector.  This  measurable  angle  is  the  angle  which  the  flow 
would  follow  if  there  were  no  viscous  effects  and  the  flow  behaved  as  a 
perfect  fluid.  Unfortunately,  a  real  flow  will  not  follow  this  ideal 
angle;  therefore,  corrections  in  the  form  of  deviation  angles  must  be 
computed.  The  accurate  determination  of  these  outlet  deviation  angles 
allows  this  analysis  to  model  a  real  flow  with  viscous  as  well  as 
secondary  flow  effects. 

The  procedure  for  calculating  the  final  deviation  angles  involves 
running  the  streamline  curvature  (SLC)  program  a  total  of  three  times. 

The  entire  procedure  is  outlined  in  the  flow  diagram  in  Figure  4.  Each 
"run"  of  the  SLC  program  adds  subsequent  terms  to  the  total  deviation. 

For  the  initial  run  of  the  program,  the  section  outlet  angles  are 
measured  and  added  to  a  deviation  term  (6..)  developed  by  Howell^. 

This  deviation  angle  is  defined  for  each  blade  section  as 


0.23  (ec)(s/c)1/2(2a/c)2 

1.0  -  (0.1/50) (0  )(s/c)1/2 
c 


(ID 
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where,  0c  is  the  section  turning  angle,  s/c  is  the  space-to-chord  ratio, 

and  a/c  is  the  distance  from  the  section  leading  edge  to  the  point  of 

maximum  camber  divided  by  the  section  chord  length. 

Once  a  converged  solution  is  obtained  for  the  flow  field  using 

Howell's  deviation,  the  axial  velocity  distributions  can  be  used  to 

calculate  an  improved  deviation  term.  The  effect  of  flow  acceleration, 

blade  camber  and  blade  thickness  can  now  be  calculated  separately.  The 

deviation  term  due  to  axial  accelerations  through  the  propeller  is  given 
7 

by  Lakshminarayana  as, 

* 

„ AVR-1.0 
A6 

TrK(c/s)(G/c+a/4)cos[(ei+62)/2][(AVR+l)2  -  4] 

AVR-1.0 

‘  8+ttK(c/s)  (a/4+G/c)cos[  (3-L+62)/2]  [  (AVR+DtanB^tanB^ 
27rK(c/s)tanB^ 

AVR  cos[(B1+82)/2]sec‘1B1  j 

8+rrK(c/s)  (a/4+G/c)cos [  (B1+B2)/2]  [  (AVR+l)tan82+2tan8^1 


+  tan82 


(12) 


where  AVR  is  the  axial  velocity  ratio  (Va2/V^) ,  K  is  the  cascade  influence 

g 

coefficient  ,  G/c  is  the  distance  from  the  chordlinc  to  the  maximum  camber 

point  of  the  section  divided  by  the  chord  length,  B^  and  82  are  the  inlet 
and  outlet  flow  angles,  respectively,  and  a  is  defined  as  the  difference 
between  the  inlet  flow  angle  (8^)  and  the  section  stagger  angle  (X). 
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Howell's  deviation  formula  considers  only  thin  blade  sections; 
therefore,  a  separate  term  due  to  thickness  effects  must  also  be  calculated. 
This  deviation  term  is  given  by  empirical  data  which  has  been  collected 

Q 

by  the  National  Aeronautics  and  Space  Administration  .  This  term  is 
calculated  from  the  equation, 

4S*  ’  •  <Vt  '  <SPlO  '  <13> 

The  terms  on  the  right-hand  side  of  Equation  (13)  are  determined  from 
graphs  of  experimental  data  based  on  section  solidity  and  inlet  flow 
angles.  Also  given  in  NASA's  development  of  deviation  terms  is  a  better 
approximation  of  the  effect  of  camber  based  on  empirical  data.  This  term 
replaces  Howell's  deviation  term.  From  graphs  of  section  solidity  and 
inlet  flow  angles,  a  value  for  the  slope,  m,  is  obtained  and  used  to 
calculate  this  deviation  term 


6  =  m  •  9 

o  c 


(14) 


Knowing  these  terms  of  the  deviation  angle  due  to  axial  acceleration 
effects,  blade  thickness  and  the  improved  camber  deviation,  one 
calculates  a  new  outlet  flow  angle  profile  from 


B 


* 

2 


B,  -  A6  +  A6  +  S 

2  o 


(15) 


This  new  outlet  flow  angle  profile  is  then  used  in  the  SLC  program  for  the 
second  run  of  this  three  part  calculation. 
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The  results  of  this  converged  solution  are  then  used  to  solve 
secondary  vorticity  equations  and  to  determine  a  deviation  term 
(A<$s)  which  is  due  to  secondary  flow  effects  as  described  by  Billet^. 
The  final  outlet  flow  angle  profile  is  obtained  by  adding  this  secondary 
flow  term  to  the  deviation  terms  thus  far  calculated  to  obtain 


8**  =  6-  -  A6'  +  A<$*  +  6  +  A<S 
2  2  os 


(16' 


This  final  outlet  flow  angle  distribution  is  then  used  as  the  second 
boundary  condition  in  the  last  run  of  the  SLC  program.  Obtaining  the 
deviation  angle  by  this  complex  procedure  enables  the  modeling  of  a  real 
fluid  flowing  through  a  real  propeller. 


Comparison  of  Results 


To  verify  the  accuracy  of  this  direct  analysis  method,  two  test 


cases  were  calculated.  In  the  first  case,  the  flow  field  was  calculated 
for  an  open  propeller  operating  at  its  design  conditions.  This  implies 
a  specific  propeller  flow  coefficient  given  as. 


(17) 


where  the  design  thrust  and  torque  are  achieved.  For  the  second  test 
case,  the  same  propeller  was  operating  at  a  10%  reduction  of  flow 
coefficient. 

These  two  cases  were  chosen  because  experimental  data  wore  available 
from  other  investigations  (References  11  and  12).  The  inlet  and  exit 
velocity  profiles  had  been  measured  in  tests  conducted  in  the  48-inch  water 
tunnel  and  the  48-inch  wind  tunnel  at  the  Garfield  Thomas  Water  Tunnel 
at  The  Pennsylvania  State  University. 
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The  comparison  between  experimental  data  and  analytical  predictions 
are  shown  In  Figures  5  through  8.  For  each  test  case  two  figures  are 
shown.  The  first  figure  In  each  case  shows  the  axial  velocity  profile 
at  the  inlet  station  to  the  propeller.  These  two  figures  show  somewhat 
similar  trends  even  though  one  is  the  design  case.  Figure  5,  and  one  is 
the  10%  off-design  case.  Figure  7.  In  these  figures  the  predicted  data 
agree  very  well  with  values  obtained  experimentally.  Some  deviation 
is  noticed  in  the  region  of  the  propeller  hub.  This  overprediction 
close  to  the  hub  is  due  to  the  inability  of  the  SLC  program  to  handle 
very  small  velocities  in  this  region.  Low  velocities  produce  program 
instabilities  and  have  to  be  handled  very  carefully.  This  stability 
problem  is  of  importance  only  on  the  first  two  or  three  streamlines  and, 
as  the  figures  verify,  has  no  effect  on  the  outer  positions  of  the  flow. 

The  second  figure  in  each  of  the  test  cases.  Figures  6  and  8,  shows 
the  velocity  components  as  a  function  of  radial  distance  from  the  propeller 
hub  at  the  propeller  exit  station.  The  design  test  case,  Figure  6,  shows 
good  correlation  between  the  predicted  and  the  experimental  axial 
velocity  data.  The  bulge  in  the  axial  velocity  profile  is  due  to  the 
contraction  of  the  streamlines  as  they  pass  through  the  propeller.  This 
contraction  is  a  result  of  the  work  done  by  the  propeller  in  the  region 
of  high  blade  loading. 

The  remaining  plot  in  Figure  6  is  a  comparison  between  the  predicted 
and  the  experimental  tangential  velocity  distribution  at  the  propeller 
exit  station.  In  this  case  the  two  curves  are  not  as  close  as  the  axial 
velocity  curves  especially  near  the  hub.  The  deviation  near  the  hub  is 
attributed  to  the  overprediction.  Figure  5,  of  the  inlet  velocity  profile  in 
this  same  region.  It  should  also  be  pointed  out  that  in  the  region 
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of  the  propeller  tip  the  flow  is  highly  three-dimensional  and  experimental 
data  may  be  expected  to  deviate  slightly  from  data  predicted  by  the 
analytical  technique. 

The  10%  off-design  case,  Figure  8,  shows  similar  trends  as  discussed 
for  the  design  case  even  though  the  absolute  magnitudes  of  the  data  are 
different.  Again,  the  correlation  between  experimental  and  analytical 
data  is  very  good  but  with  small  deviations  near  the  hub  of  the  tangential 
velocity  curve. 


IV.  SUMMARY 

A  numerical  solution  to  the  direct  problem  has  been  developed  for 
calculating  the  flow  field  of  an  open  propeller  using  the  streamline 
curvature  method.  Results  were  obtained  which  closely  model  the  flow 
by  implementing  the  streamline  curvature  technique  to  solve  the  equations 
of  motion  through  the  turbomachine.  By  making  adjustments  to  the  outlet 
flow  angles  to  allow  for  viscous  and  secondary  flow  effects,  good 
correlation  between  experimental  measurements  and  theoretical  predictions 
have  resulted.  The  degree  of  accuracy  which  has  been  experienced  thus 
far  indicates  that  this  method  can  be  a  valuable  tool  in  the  performance 
predications  of  a  propeller  operating  at  other  than  design  conditions. 
Consequently,  this  method  has  the  capacity  to  give  the  engineer  a  clear 
picture  of  how  a  particular  propeller  will  act  at  off-design  conditions 
without  going  through  the  time  and  expense  of  experimental  testing. 
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BOUNDING  STREAMSURFACE 


Figure  la  -  Velocity  Components  and  Coordinate  System  as  Seen  in  a  Side 
View  of  Flow  Field  (Meridional  Plane) 
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Figure  lb  -  Velocity  Components  as  Seen  in  a  Top  View  of  Flow  Field 


Diagram  of  SLC  Computational  Procedure 
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Figure  3  -  Propeller  Section  Geometry 
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Figure  4 


A  Flow  Diagram  of  the  Procedure  Used  to  Calculate  the  Deviation 
Angles 
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Figure  6 


-  Comparison  of  Measured  and  Predicted  Velocity  Profiles  at 
Propeller  Exit  Station  for  Design  Conditions 
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Figure  7 
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-  Comparison  of  Measured  and  Predicted  Velocity  Profile  at 
Propeller  Inlet  Station  for  10%  Off-Design  Conditions 
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Figure  8  -  Comparison  of  Measured  and  Predicted  Velocity  Profiles  at 
Propeller  Exit  Station  for  10%  Off-Design  Conditions 
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